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Abstract 


Nomenclature 


The accuracy of a nonintrusive high-angle-of-attack 
flush airdata sensing (HI-FADS) system was verified 
for quasi-steady flight conditions up to 55° angle of 
attack during the F- 18 High Alpha Research Vehicle 
(HARV) Program. The system is a matrix of nine 
pressure ports arranged in annular rings on the air- 
craft nose. The complete airdata set is estimated us- 
ing nonlinear regression. Satisfactory frequency re- 
sponse was verified to the system Nyquist frequency 
(12.5 Hz). The effects of acoustical distortions within 
the individual pressure sensors of the nonintrusive 
pressure matrix on overall system performance are 
addressed. To quantify these effects, a frequency- 
response model describing the dynamics of acoustical 
distortion is developed and simple design criteria are 
derived. The model adjusts measured HI-FADS pres- 
sure data for the acoustical distortion and quantifies the 
effects of internal sensor geometries on system perfor- 
mance. Analysis results indicate that sensor frequency 
response characteristics vary greatly with altitude, thus 
it’s difficult to select satisfactory sensor geometry for 
all altitudes. The solution uses presample filtering to 
eliminate resonance effects, and short pneumatic tub- 
ing sections to reduce lag effects. Without presam- 
ple signal conditioning the system designer must use 
the pneumatic transmission line to attenuate the reso- 
nances and accept the resulting altitude variability. 
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A c 

cross-sectional area of pneumatic 
tubing 

A(a) 

intermediate coefficient of integration 

B(s) 

intermediate coefficient of integration 

BVP 

boundary value problem 

C p 

pressure coefficient 

c 

sonic velocity 

D 

pressure tubing diameter, in. 

ESP 

electronically scanned pressure 

ESP101 

OP clock angle, 0° cone angle, 
HI-FADS pressure 

e 

complex exponential 

T 

Fourier transform operator 

F(. . .) 

nonlinear functional (potential flow) 

FADS 

flush airdata sensing 

FFT 

fast Fourier transform 

FRP 

fiberglass reinforced plastic 

f 

input frequency 

HARV 

high alpha research vehicle 

HI-FADS 

high-angle-of-attack flush airdata 
system sensing 

Hp 

pressure altitude 

i 

HI-FADS pressure port index 


complex constant, (— l) 1 / 2 

L 

length of pneumatic tubing, ft 

C 

Laplace transform 

Ml 

Mach number from output pressures 

Mo 

Mach number from input pressures 



m 

summation index 

N 

number of HI-FADS pressure ports 

P 

ambient pressure 

Pl(D 

Fourier domain 

Pl(s) 

pressure function output (frequency 
domain) 

Po(/) 

Fourier domain 

Po(s) 

pressure function input (frequency 
domain) 


second order transfer function 

P(z,s) 

general pressure function (frequency 
domain) 

PCM 

pulse code modulation 

PSD 

power spectral density 

Poo 

aircraft static pressure 

Pi 

pressure at/th HI-FADS orifice 

PL(t) 

pressure function output (time domain) 

po(0 

pressure function input (time domain) 

p(z,0 

general pressure function (time 
domain) 

p(&,^») 

pressure value at FADS coordinates, 
4>i> 

Q 

incompressible dynamic pressure 

Qc 

compressible dynamic pressure 

R 

acoustic resistance of pneumatic tubing 

s 

Laplace operator 

T(x, t) 

temperature function 

t 

time coordinate 

U(x,t) 

velocity function 

V 

enclosed transducer volume, in 3 

V 

frequency domain transfer function 

V 

time domain describing function 

X 

spatial coordinate 

a 

angle of attack (def. 1) 

a 

complex plane frequency variable 
(def. 2) 

(*L 

angle of atttack from output pressures 

ao 

angle of atttack from input pressures 

P 

angle of sideslip 

6 

coefficient of higher order harmonics 


£ 

HI-FADS calibration parameter 

c 

damping ratio 

X 

HI-FADS cone angle coordinate 

p 

dynamic viscosity 

po 

inital or ambient density 

T 

dummy variable of integration 
(def. 1) 

T 

acoustical time lag (def. 2) 


HI-FADS clock angle coordinate 

Wn 

natural frequency 


Introduction 

Contemporary design requirements for some high- 
angle-of-attack and most hypersonic or low observable 
aircraft do not allow for airdata measurement using 
conventional noseboom instrumentation. Instead, re- 
quirements for these classes of vehicles mandate that 
airdata be measured nonintrusively, that is, without 
placing a sensing probe directly into the flow stream. 
One very promising class of nonoptical, nonintru- 
sive airdata systems, the flush airdata sensing (FADS) 
system 1-3 uses a matrix of surface pressure measure- 
ments to infer free-stream airdata indirectly. 

A particular system designed for high-angle-of- 
attack measurements, the high-angle-of-attack flush 
airdata sensing (HI-FADS) system, Fig. 1, was veri- 
fied to be accurate for quasi-steady conditions up to 
55° angle of attack during phase one of the F-18 high 
alpha research vehicle (HARV) flight test program. 
Preliminary test results were reported in Ref. 1. The 
HI-FADS system consists of a matrix of pressure ports 
arranged in an nular rings on the aircraft nose, and es- 
timates the complete airdata set using potential flow 
modeling and nonlinear regression. 

The system performed well for steady or moderate 
maneuvering flight conditions and computations were 
performed at a rate of 25 samples/sec. Satisfactory fre- 
quency response was verified up to the system Nyquist 
frequency of 12.5 Hz. 1 For rapid maneuvering and un- 
steady flight regimes, where higher frequency airdata 
could be utilized, there is some question as to how 
well the HI-FADS system would perform. There are 
two primary reasons for this uncertainty: (1) unmod- 
eled, unsteady aerodynamic effects, and (2) acousti- 
cal distortion within the individual pressure sensors of 
the HI-FADS surface pressure matrix. Research pre- 
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sented in Ref. 4 addresses the unmodeled, unsteady 
aerodynamics. This paper will address the effects of 
acoustical distortion. 

Background information on the physics of pressure 
sensor acoustical distortion is presented first Next 
an accurate model which describes the dynamics of 
acoustical distortion will be presented. The model, 
which has been verified extensively 5-8 will be briefly 
verified for the HI-FADS configuration by compar- 
ison to laboratory test results. Then system design 
criteria will be developed from the model. Applica- 
tion of these criteria to the HI-FADS system design 
will be discussed. Using the techniques of spectral 
deconvolution 7 ~ 9 the model is used to adjust mea- 
sured HI-FADS pressure data for the effects of acous- 
tical distortion and the results are merged with other 
high-frequency surface pressure flight data to synthe- 
size a representative surface pressure data set. Using 
the synthesized surface pressure data set for sensor in- 
puts, the effects of various internal sensor geometries 
will be analyzed and qualitative comparisons will be 
presented. Various system design difficulties will be 
discussed. 

Airdata Measurement System 

The HI-FADS configuration has a simple hardware 
arrangement with the basic fixture being a fibeiglass- 
reinforced-plastic (FRP) cap, mounted on the nose of 
the F-18 HARV. A set of 25 0.06-ia inside diame- 
ter pressure orifices, arranged in 4 annular rings, were 
drilled in the nosecap. Flight tests conducted during 
phase one of the HARV Program, indicated that air- 
data could be measured satisfactorily using a subset of 
nine pressure ports from two of the rings. All results 
presented in this paper will use the 9-port configura- 
tion. Pressures at the nosecap were sensed by a multi- 
ple transducer, electronically scanned pressure (ESP) 
module remotely mounted on a structural bulkhead 
in the aircraft nose cavity. The ESP module, which 
consists of differential transducers, was referenced to 
a single, high-accuracy absolute pressure transducer 
mounted in the aircraft nose cavity. High-frequency 
dynamics in the reference pressure were damped by 
a 50.0 in 3 reference tank also mounted in the aircraft 
nose. The temperature environments of both the ESP 
module and the reference transducer were controlled 
by wrapping the units in heater blankets to maintain an 
operating temperature of approximately 1 1 0 °F. For the 
ESP module the estimated accuracy is ±1.0 lb/ft 2 and 


for the reference transducer, the estimated accuracy is 
±0.50 lb/ft 2 (Ref. 1). The outputs from both the ESP 
module and the reference transducer were digitized by 
pulse code modulation (PCM) and telemetered to the 
ground at 25 samples/sec. Because of telemetry limi- 
tations, higher sample rates were not practical for the 
HARV flight tests. 

Since the ESP module was electronically scanned 
(multiplexed), any information with significant power 
which lies beyond the Nyquist frequency of the scan 
(12.5 Hz) will be aliased to a lower frequency and 
will distort the measured pressure values. Since the 
transducer output is multiplexed — the equivalent of 
time sampling — posttransduction signal conditioning 
cannot be used without eliminating information in the 
band of interest. For the HARV HI-FADS installation, 
the sensing ESP pressure transducer was located re- 
motely from the surface, and 8-ft lengths of 0.06-in. - 
diameter flexible pneumatic tubing were used to trans- 
port pressure from the surface to the transducer. Each 
ESP transducer has an internal volume of less than 
0.01 in 3 . As discussed in Ref. 1, the previously men- 
tioned tubing geometries were selected to give the sen- 
sors low-pass characterisitcs and thus serve as crude 
antialiasing filters. 

For the remainder of the paper, the term “trans- 
ducer” will be used to describe the actual ESP physical 
transduction device. The term “sensor” will be used 
to describe the transducer and all associated pneumatic 
tubing, pressure fittings, and so on. The term “system” 
will be used to describe the collection of all nine HI- 
FADS measurement sensors, the pressure reference 
system, and the associated data acquisition, teleme- 
try, and storage hardware. A schematic of the HI- 
FADS hardware arrangement is presented in Fig. 1. A 
schematic of the internal geometry of the HI-FADS 
pressure sensing system is presented in Fig. 2. More 
detail concerning the HI-FADS system hardware and 
the HARV research data acquisition system may be 
found in Ref. 1. 

The measured pressure data, p( fa, A<) , were related 
to the free-stream airdata quantities using a modified 
potential flow model 1 

p( ^i> ^i) = F*( Of, 0, q c , Poo i ^i,^i,£),r— 1, N 

and all nine pressure observations were used simulta- 
neously to estimate the four airdata parameters: com- 
pressible dynamic pressure (q c ), angle of attack (a). 
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angle of sideslip (/?), and static pressure (p*), using a 
least-squares criterion and nonlinear regression. Here 
<f>i and \i are the clock and cone angles of the ith 
HI-FADS pressure port, e is a calibration parameter, 
and N is die number of ports (nine in this case). Using 
these four basic airdata parameters, most airdata quan- 
tities of interest may be calculated directly. The use 
of an overdetermined (more observations than states) 
analysis minimizes the effects of errors in any single 
pressure measurement. The resulting algorithm, here- 
after referred to as the “HI— FADS algorithm,” is ro- 
bust and accurate. Excellent steady-state system per- 
formance has been verified over the entire subsonic 
Mach number range and up to 55° angle of attack. Spe- 
cific details concerning the algorithm implementation 
may be found in Ref. 1 . 

Sensor Acoustical Model 

Referring to Fig. 2, the sensor configuration is mod- 
eled as a straight cylindrical tube with an axisymmet- 
ric volume, which models the internal volume of the 
pressure transducer attached to its downstream end. 
The tube is considered to be of constant diameter ( D ) 
and length (L) with a longitudinal coordinate (z) de- 
fined from the upstream -end of the tube. Pressure, den- 
sity, temperature, and flow velocity within the tube are 
modeled as one-dimensional functions of longitudinal 
position and time, and are represented by the symbols, 
p(z,f), T(x,t), and U(x,t). For this analysis the 
flow velocity, U(x,t), although actually a function of 
radial position within the tube, is approximated by its 
radial average. The input pressure, the value at x = 
0, is represented by the symbol po(t). The measured 
pressure at the transducer, the value at x = L, is repre- 
sented by the symbol pi( t ) . 

Pressure variations at the surface propagate as lon- 
gitudinal waves through the connective tubing to the 
transducer. The wave propagation is damped by fric- 
tional attenuation along the walls of the tubing. When 
the wave reaches the downstream -end of the tubing, it 
reflects back up the tube and may either damp or am- 
plify incoming pressure waves. The frictional damp- 
ing and wave interference produce a magnitude distor- 
tion and phase delay in the measured pressure signal. 
The magnitude of the distortion and phase delay are 
dependent on the internal geometry, surface crossflow 
on the sensor pressure port, and the free-stream flight 
conditions — primarily altitude. The speed of the lon- 
gitudinal wave is somewhat slower than the local sonic 


velocity but considerably faster than the internal flow 
velocity. Analytical and experimental results 5-9 indi- 
cate that for modem, low-volume installations, inter- 
nal flow velocities are considered to be small. 

The dynamics of the acoustical wave propagation 
are described accurately by a boundary value prob- 
lem (BVP) derived from the Navier-Stokes equations 
of momentum and continuity. 5,9 The acoustical wave 
propagates according to 


d 2 p(x,t) R dp(x,t) _ 2 c^pOM) 
dt 2 + po dt C dx 2 

Subject to the boundary conditions 


d 2 p(L,t) , 
dt 2 

+ 


R dptj^t) 
po dt 


A c < c 2 

V 


{ 


dp(s,t) I 
dx J X=L 


= 0 


(lb) 


p(0 ,t) = po(t) ( 1c) 

and the initial condition 

p(x,0) = P a (Id) 

In equations (la) and (lb), R is the acoustical resis- 
tance of the system and as developed in Ref. 7 for lam- 
inar flow 



and A c is the cross-sectional area of the tube. By evalu- 
ating the acoustical resistance based on Blasius friction 
law 10 the BVP may be easily extended to describe tur- 
bulent flow. This extension is derived in detail in Refs. 
5 and 8. In general, for the case of arbitrary pressure 
inputs the BVP cannot be solved analytically. How- 
ever, numerical techniques for solving the BVP for ar- 
bitrary inputs are developed in Refs. 5 and 8. The BVP 
can be solved in closed form, however, for the case of 
sinusoidal inputs — commonly referred to as the “fre- 
quency response” of the model. 

Frequency Response of the Acoustical Model 

The frequency response solution is developed by 
first taking the Laplace transform of equations ( 1 a) and 
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(lc) and removing the constant initial conditions by Letting a — for the Fourier transform and eval- 

subtracting the indal pressure value. The result is uating at x = L, simplifying and reducing gives, 


n , s 2 d 2 P(x,s) 

P(x,s) = C ( 2a ) 


(* 2+ po _ dx* 

Age 2 f dPjXjj) 


2 fl \ n/ , . A e c* f dP(x,a ) 1 _ A 

* + — a Pl(s) + —rr- \ jr \ = 0 

\ Po / V V ox ) x= i 

(2b) 

Po(s) = £[po(t)l 
Integrating with respect to x gives 

P(x,s) = A(s)t^ /c + B(s)e~^ x/c (3a) 


and differentiating equation (3a), 


dP(x,s) 

dx 



e \fcx/c _ B(g) e -VZ*/c 


(3b) 


where. 


2 ^ 

a = s + — 3 
P 0 


Evaluating P(x,s ) (eq. (3a)) at x = 0, and 

(eq. (3b)) at x = L and solving for A(s), and B(s), 

then 


A(s) = 


yPo(s) [l — ] e-v^/ c 

cosh [ y/aL/c] + sinh [ y/aL/ c ] 


(3c) 


and 


B(s) = P 0 (s) 


^Po(a) [l ~ n§] e~' /SL,e 

cosh [y/aL/c] + sinh [ ^/aL/c] 

(3d) 


Substituting these expressions into equation (3a), 


PlU) 1 

Po(f) cosh [ y/aL/c] + sinh [y/aL/c] 

(4) 

where, 

a= -(4ir 2 / 2 ) + 

Equation (4) is a simple frequency response model in 
which the ratio of the complex spectra is given explic- 
itly as a function of the sensor geometry, the initial flow 
density (usually taken as ambient), and the frequency 
of the input sinusoid. Results similar to equations (3e) 
and (4) are derived in Refs. 11 and 12. 

Verification of the Acoustical Model 

The acoustical model of equations (la)-(ld), has 
been extensively verified by comparisons of numeri- 
cal solutions to both laboratory and flight data. The 
results of this verification are presented in detail in 
Refs. 5-8. The closed form solution, equation (4), has 
been verified by extensive comparisons to both numer- 
ical solutions and laboratory data. Figure 3 presents 
a sample frequency response comparison. Depicted 
in Fig. 3 are frequency response comparisons of lab- 
oratory data and values generated using equation (4). 
In this figure, data were obtained using the HI-FADS 
pressure sensor geometry (L = 8 ft, D = 0.06 in., and 
V — 0.01 in 3 ) and an ESP transducer. The data pre- 
sented were obtained at 2300- ft altitude. Although the 
lab data indicate that the acoustical model is slightly 
less damped than the actual sensor, considering that the 
analyses of equations (la)-(ld) do not consider tubing 
roughness or constrictions, the comparisons are quite 
favorable. Detailed descriptions of the laboratory and 
flight experimental setups may be found in Refs. 5-8. 


n_.) f 

[ cosh [ y/aL/c ] + sinh [ y/aL/c ] 

x i r . t/. /z~] 


ft(«) 


) e \/5x/c 

_ 

iPbU) [l - ^f] «-' /SWc ' 

™ t/ /r. " _ 


coshf v/aL/cl 


(3c) 
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Convolution and Deconvolution 
of Pressure Signals 

Equation (4), describes the behavior of the pres- 
sure sensor as a function of the frequency of the input 
sinusoid. It has multiple resonances, and can be ap- 
plied to arbitrary input data by decomposing the input 
signal into its Fourier components by way of Fourier 
Transform techniques 13,14 . Mathematically, the phys- 
ical process may be modeled by the convolution 


piXt) = [ v(t-T)poCr)dT (5a) 

Jo 

where, v( t) , is a describing function for the sensor dy- 
namics. Evaluating the Fourier transform of equation 
(5a), then 


fiplU)] = t 



- r)po(r)dr 


where. 


= V(/)Po(/) 
(5b) 


V(/) = 


PlW 

Po(/) 

! (5 C ) 

cosh [^/aL/c} + ^^sinh [y/aL/c\ 


and is referred to as the sensor “transfer function.” 
The time domain equivalent of the output spectrum, 
Pl( /) , may be generated by numerically evaluating 
the inverse Fourier transform. 


The numerical transform operations are easily 
mechanized using fast Fourier transform (FFT) tech- 
niques. The specific analysis procedure, referred to 
as “spectral convolution,” is depicted schematically 
in Fig. 4. Since finite record lengths are being used, 
to prevent spectral power leakage 13-15 care must be 
taken to appropriately window the input time history 
data. For this analysis a cosine-taper window was 
used. This procedure is shown in block 1 of Fig. 4. 
The data are then transformed using the FFT pro- 
cedure (block 2). To eliminate noise introduced by 
the discrete transformation, as well as noise in the 
original data, the resulting spectra are smoothed us- 
ing frequency averaging. For this analysis, a two- 


pass recursive implementation 14 was used.To elimi- 
nate frequency-shift distortions, the first pass was per- 
formed forward in frequency and the second pass was 
performed backward in frequency (block 3). The re- 
sulting smoothed spectrum is convolved to account 
for acoustical distortion by multiplying by the trans- 
fer function of equation (5c). This procedure is shown 
in block 4. Finally the results are inverse transformed 
to give the time history (block 5). 

The procedure is easily reversible, and up to the 
Nyquist frequency of the system, it allows the input to 
be inferred from a measured output The reverse pro- 
cess is performed identically, except that the inputs and 
outputs are reversed and in block 4 the smoothed spec- 
trum is divided (instead of multiplied) by the transfer 
function of equation (5c). This process is referred to 
as “spectral deconvolution,” and is depicted schemat- 
ically in Fig. 5. 

Formulation of Qualitative Design Criteria 


If the denominator terms of equation (4) are ex- 
panded in a Taylor’s series 16 


cosh 



j 


(6a) 


and 


sinh 



2 ™ +1 a ( 2n .+ l)/2 

(2m + 1)! 


} (6b) 


and terms in like powers of a are collected, then 
equation (4) may be expressed as an infinite series of 
harmonics 


PiAs) 

Po(s) 


1 


l + 


1 

on 

+ V 

[(*ri 

i 

l 

(2m)! 

AcC 

(2m— 1) ! 

} 


a 1 ' 


(6c) 


evaluating 

2 R 

a = s + — s 
P o 
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and regrouping terms then 


of the full-wave model. Thus, if the resonance result- 
ing from the dominant harmonic is eliminated by ad- 


PUs) 

Po(s) 


1 







(6d) 


If all terms (in s) of equation (6d) of an order greater 
than 2 are neglected, then a reduced-order model 
which describes the behavior the primary or dominant 
harmonic results, 

(PlC s)\ 

VPo(s) / \h 


herence to the design rules, either by geometrical de- 
sign or presample low-pass filtering, then the sensor 
will not resonate. 



Equation (6e) is a second-order filter of the form 


( Pl(s) \ _ l 

\Po(s)J lh -W + 2-fa+l 


(6f) 


with a natural frequency given by (assuming laminar 
flow) 




m 


{[*] [> + i&M'+S]} 

a damping ratio given by 

„ 


(6g) 


c + 


2 po 


ill LV 

2 c 2 A-c 2 


(6h) 


and a steady-state time constant given by 


Equation (6h) indicates that the sensor damping ra- 
tio is inversely proportional to flow density. Since den- 
sity varies inversely with altitude, a sensor which is un- 
derdamped at low altitude may be highly overdamped 
at high altitudes. Again, use of the design criteria to 
describe the first harmonic allows for the geometry of 
the pressure sensor to be quickly tailored to a particular 
frequency band or operating altitude. 

Similar numerical evaluations of equation (4) indi- 
cate that the time lag induced by the first harmonic ac- 
counts for most of the pressure sensor measurement lag 
and that equation (6i) is a good indicator of the over- 
all sensor time lag. Thus, based on the behavior of the 
first harmonic, the design criteria allow the general dy- 
namical behavior of the pressure sensor to be quickly 
and easily predicted. 

Results and Discussion 


2 c 128 fiL \v + 

T -g- 1 (6i) 

The time constant of equation (6i), is identical to 
the value predicted by the first-order analysis of 
Ref. 17 — long considered to be the best available pre- 
diction of acoustically-induced pressure measurement 
lag. Taken together, equations (6g) through (6i) con- 
stitute an important set of qualitative design criteria. 

Extensive lab tests 6 and numerical evaluation equa- 
tion (4) (full-wave model) indicate that peak magni- 
tudes of the higher order harmonics are smaller than 
the peak magnitude of the dominant harmonic. Fur- 
thermore, evaluation of equation (6g) (reduced-order 
model) for many geometries indicates that the natu- 
ral frequency and peak magnitude of the second-order 
model are a good approximation of the first harmonic 


The concepts discussed previously will now be illus- 
trated. First, the frequency response of several sensor 
geometries and altitudes will be illustrated using the 
exact solution of equation (4) and the design rules of 
equations (6g)-(6i). Next, the time response will be 
illustrated using data synthesized from HI-FADS and 
other flight data. Power spectra, time history compar- 
isons, and correlation analyses will be presented. 

Application of Design Criteria 

Figure 6 presents the calculated frequency response 
of a typical HI-FADS pressure sensor at 20,000- ft alti- 
tude. Here the frequency response is flat to 10 Hz and 
has an attenuation of approximately 1 dB at 25 Hz. 
Beyond this frequency, the magnitude rolls-off more 
steeply. Note that the 2 small resonance harmonics 
which occur at approximately 90 and 1 80 Hz have been 
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significantly attenuated by the pneumatic tubing. Be- 
cause the harmonics are attenuated, their resonances 
will not diston the measured pressure signal in a sig- 
nificant manner. As calculated using equation (6h), the 
approximate time delay of the pressure sensing system 
is 15 msec. The HI-FADS system was designed for an 
operating altitude range between 10,000 and 40,000 ft, 
and the internal pressure sensor geometries were 
selected by adherence to the design rules of equa- 
tions (6gM6i) as a compromise between these two al- 
titude extremes. At 10,000-ft altitude the pressure sen- 
sors were slightly underdamped, while at 40,000-ft al- 
titude the sensors were somewhat overdamped. 

In figure 6, both of the harmonics lie beyond the 
12.5 Hz Nyquist frequency of the HI-FADS system. 
Recall that the sample rate was 25 Hz. If the har- 
monics were not attenuated, then the resonated noise 
would have been aliased to the response band of 
the HI-FADS system. In other words, using signifi- 
cantly wider or shorter lengths of tubing at this alti- 
tude would have resulted in degraded pressure mea- 
surements. This effect is depicted in Fig. 7, where, 
instead of using 8-ft tubing lengths, the frequency re- 
sponse resulting from use of 1-ft tubing sections is 
evaluated. Here a very strong harmonic is present at 
approximately 215 Hz. In designing a HI-FADS type 
of airdata system, in which acoustical resonance may 
be present, great care must be taken to range the pres- 
sure sensors to accommodate any magnitude amplifi- 
cation which may be induced by resonance. This is 
especially true if fine-scale differential measurement 
transducers are used. If a resonance condition devel- 
ops and causes noise to be so greatly amplified that the 
transducer saturates, then this nonlinear process cannot 
be accounted for and the data will be nonrecoverable. 

If significantly thinner or longer lengths of tubing 
were used, or if the system was operated at a sig- 
nificantly higher altitude, then the response of the 
measurement sensors would have been greatly atten- 
uated and unacceptable phase delays would have re- 
sulted. This effect is illustrated in Fig. 8, where the 
response of an L = 8 ft, D = 0.06 in. sensor is eval- 
uated at 65,000-ft altitude. In this case the system is 
highly overdamped and the approximate sensor time 
lag, evaluated from equation (6h), is 120 msec. Such 
phase delays are unacceptable for dynamic use. 

As a final illustration, if 1-ft tubing lengths are used 
at 65,000-ft altitude, then the system frequency re- 


sponse again is seen as satisfactory. This case is illus- 
trated in Fig. 9. Here, the frequency response of the 
pressure measurement sensor is flat to beyond 35 Hz 
and no resonant harmonics appear. The approximate 
time lag is 2 msec. 

Because the sensor frequency response characteris- 
tics vary so greatly with altitude, it is difficult to se- 
lect any single sensor geometry which gives acceptable 
acoustical performance at all altitudes. Fortunately, 
the problem can be overcome by use of appropriate 
presample signal conditioning. Using the design rules 
of equations (6h) and (6i) to identify the location of 
the primary harmonic, analog filters may be designed 
to attenuate the resonances. At the same time the fre- 
quency response characteristics of the lower frequency 
band will be preserved. If these filters are coupled to 
transducers mounted with short sections of pneumatic 
tubing, then the result is end-to-end pressure sensors 
which eliminate the resonance effects while leaving 
the lower frequency band data unaltered. Again, cau- 
tion must be used in scaling the transducer range to 
accommodate resonance. If the resonance causes the 
transducer to saturate, then when the transducer output 
is run through the low-pass filter the nonlinear distor- 
tion caused by the saturation process will result in a 
gain offset and the measured signal will be biased. 

As an example, assume a second-order, low-pass fil- 
ter with a natural frequency of 50 Hz and a damping 
ratio of 0.7071 (Butterworth filter) is coupled with the 
L = 1 ft, D = 0.06 in., V - 0.01 in 3 sensor geometry at 
20,000-ft altitude. For this case the frequency response 
would appear as in Fig. 10. Comparing this frequency 
response to Fig. 7, notice that although the resonance 
peak at 215 Hz has been significantly attenuated, the 
frequency response of the system below 40 Hz is flat. 
The coupled Butterworth filter adds only 4.5 msec of 
steady-state time delay to the measured pressure sig- 
nal. Consequently, in the 40 Hz and below frequency 
band, the pressure measurements are essentially free 
of acoustical distortion. 

The frequency response of the same tubing-low- 
pass filter configuration at 65,000-ft altitude is pre- 
sented in Fig. 11. Here again, up through 40 Hz the 
frequency response curve is flat with an attenuation of 
only 1.5 dB at 40 Hz. Pressure inputs in this band 
would be transmitted without significant acoustical 
distortion. Comparing Fig. 1 1 to Fig. 9, in the 40 Hz 
and below region the magnitude curves show no ap- 
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preciable difference and the additional phase lag intro- 
duced by the low-pass filter — less than 60° — is negli- 
gible. The time delay caused by the low-pass filter at 
40 Hz is 4.2 msec. 

At this point, a word of caution is appropriate 
concerning the use of multiplexed or electronically 
scanned pressure transducers. If multiplexed sensors 
such as the ESP module are used, then the scan rate 
must be high enough to clearly capture any significant 
tubing resonances which occur. Recall that multiplex- 
ing is the equivalent of time-sampling, and any har- 
monic which lies beyond the Nyquist frequency of the 
scan will appear at a much lower frequency. Conse- 
quently, postmultiplexing signal conditioning, either 
digital or analog, cannot attenuate the aliased reso- 
nances. If high scan rates are not possible, or pre- 
sample signal conditioning is unavailable, the system 
designer must use the pneumatic transmission line to 
attenuate the resonances and accept whatever altitude 
variability occurs. For a HI-FADS type of airdata sys- 
tem, if a wide altitude range is required and presam- 
pling filtering cannot be performed, then redundant ex- 
ternal ports with internal geometries tailored to specific 
altitudes can be used. 

Flight Data Evaluation 

The spectral distortion effects described in the pre- 
vious section will now be illustrated using flight-based 
synthesized data. To accurately simulate the actual 
presample inputs to the HI-FADS pressure measure- 
ment system, a high sample rate (500 samples/sec) sig- 
nal was synthesized from combinations of HI-FADS 
and other flight data. First, measured HI-FADS pres- 
sure data from a pushover-pullup (POPU) maneuver 
obtained at a Mach number of approximately 0.25, an 
angle of attack varying from 10° to 45°, and an al- 
titude of 20,000 ft, were deconvolved using the pro- 
cedure outlined in Fig. 5, assuming a tubing length of 
8 ft, a diameter of 0.06 in. and a transducer volume of 
0.01 in 3 . The deconvolution was performed out to the 
Nyquist frequency of the system, 12.5 Hz. 

The deconvolved HI-FADS pressure data were in- 
herently upper band-limited at 12.5 Hz. The data 
were merged with 500 samples/sec surface pressure 
flight data (obtained under similar conditions) by a 
piezoelectric microphone mounted flush to the air- 
craft surface. 6,8 This process was performed to re- 
construct the high-frequency end of the pressure spec- 
tra. For this analysis, the surface pressure data were 


lower band-limited using a high-pass filter with a break 
frequency of 12.5 Hz. The resulting time histories 
were spliced using a complementary filter in which 
the transforms of both the deconvolved HI-FADS data 
and the high-pass-filtered surface pressure data were 
evaluated. The spectra were scaled and added, and the 
resulting spectrum was inverse transformed to give a 
500 samples/sec time history. The merging process, 
essentially a complementary filter, was performed for 
each of the nine HI-FADS pressure measurements. A 
flowchart of the complementary filtering process is de- 
picted in Fig. 12. 

Once the HI-FADS pressure data were deconvolved 
and merged with the surface measurements, corre- 
sponding airdata values were evaluated using the 
HI-FADS algorithm. For the remainder of this anal- 
ysis, the deconvolved and merged pressure data will 
be referred to as the “input” pressure data, and the re- 
sulting airdata values will be referred to as the “input” 
airdata. Input Mach number will be represented by the 
symbol. Mo . Input angle of attack will be represented 
by the symbol, ao . The resulting Mach number and 
angle-of-attack time histories are presented in Fig. 13. 

To illustrate the effects of various sensor geome- 
tries, the input spectra will now be convolved using the 
method outlined in Fig. 4. The resulting pressure data 
will be referred to as the “output” pressure data. The 
convolved data set will then be passed through the HI- 
FADS algorithm to generate a corresponding airdata 
set which will be referred to as the “output" airdata. 
Output Mach number will be represented by the sym- 
bol, Mi. Output angle of attack will be represented by 
the symbol, ai. 

First the input pressure data were convolved to sim- 
ulate the effects of the HI-FADS configuration with 
8 ft of 0.06-in.-diameter pneumatic tubing and a trans- 
ducer volume of 0.01 in 3 at 20,000-ft altitude. Com- 
parisons of the input and output spectra are depicted in 
Fig. 14. Depicted in Fig. 14(a) are the spectral densi- 
ties obtained from the nosetip pressure port, ESP 1 0 1 . 1 
The ESP101 data are qualitatively representative of all 
the HI-FADS pressure data. Depicted in Fig. 14(b) 
are the corresponding angle-of-attack spectral densi- 
ties. Note that output pressure data, Pi, show sig- 
nificant power losses at higher frequencies, however, 
the output angle of attack, ai, shows only a minimal 
power loss when compared to the input angle of at- 
tack, ao. This effect is a result of the HI-FADS algo- 
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rithm least-squares regression eliminating noncoher- 
ent, high-frequency noise from the set of nine input 
HI-FADS pressure measurements. The output pres- 
sure data already have been attenuated at higher fre- 
quencies by the pressure tubing (convolution) and the 
resulting input and output angle-of-attack spectra are 
neatly identical. 

Corresponding Mach number and angle-of-attack 
time histories are presented in Fig. 15. Here, the time 
scale has been expanded to illustrate time delays and 
magnitude differences. Note the time lag when com- 
paring the output signals, Ml and o*L> and the input 
values. Mo and ao- Cross-correlation analyses be- 
tween the input and convolved data values indicate that 
for this geometry the output angle-of-attack signal is 
acoustically lagged by approximately 15 msec. The 
cross-correlation data are presented in Fig. 16. 

The effect of resonance is now illustrated by con- 
volving the input pressure data, assuming a sensor tub- 
ing length of 1 ft. The diameter is assumed to remain 
at 0.06 in. and the transducer volume is assumed to re- 
main at 0.01 in 3 . Resulting spectra are presented in 
Fig. 17. Presented are the spectra for ESP101, and 
angle of attack. Note that the output pressure is res- 
onated at high frequencies. Also presented are the cor- 
responding angle-of-attack spectra. The regression of 
the HI-FADS algorithm does not remove the resonated 
noise, which is coherent, from the input pressure data. 

Corresponding Mach number and angle-of-attack 
time histories are compared in Fig. 18. Again the time 
scale has been expanded to illustrate the effects. Note 
that both the output Mach number and angle-of-attack 
values are noisy. The resonance-induced acoustical 
noise results in an angle-of-attack noise band of ap- 
proximately one-half of a degree and the correspond- 
ing Mach number noise band is 0.004. 

Since no flight data for extreme altitudes were ob- 
tained for phase one of the HARV flight tests, the ef- 
fects of very high altitudes will now be illustrated by 
extrapolating flight data to 65,000-ft altitude from the 
maneuver of Fig. 13. The input pressure data were 
extrapolated to 65,000 ft (kft = ft x 10 3 ) by assuming 
that 

(CpQ)6Skft _ ( P; - Po>65 Jfc/t _ j 
(Cp Qhokft (P; - Pahokft 


Resulting input Mach number and angle-of-attack time 
histories, generated by the HI-FADS algorithm, are 
presented in Fig. 19. 

Using the method of Fig. 4, the synthesized input 
data are convolved assuming a sensor tubing length 
of 8 ft. The diameter and transducer volume remain 
at 0.06 in. and 0.01 in 3 , respectively. The result- 
ing pressure data are passed through the HI-FADS 
algorithm to give an output airdata set. The result- 
ing angle-of-attack spectra and time histories are com- 
pared in Fig. 20. Cross-correlation data are presented 
in Fig. 21. Note that the output angle-of-attack data 
are highly lagged, with an approximate time delay of 
100 msec. This acoustical delay is unacceptable for 
dynamic usage. The 8-ft tubing section HI-FADS 
geometry, although giving acceptable performance at 
20,000 ft, induces unacceptable phase lag at high 
altitudes. 

As discussed earlier, the resonance and lag difficul- 
ties can be overcome by use of presample signal con- 
ditioning. The performance of the 1-ft tubing length 
configuration coupled with the 50 Hz Butte rworth fil- 
ter will now be analyzed. Figure 22 presents the spec- 
tral comparions for data obtained from the 20,000- 
ft-altitude maneuver. Notice that most resonance- 
induced noise has been removed from the pressure 
and angle-of-attack signals, and the input and output 
angle-of-attack spectra are nearly identical. Figure 23 
presents the corresponding time history data. Clearly, 
the resulting time history traces are nearly indistin- 
guishable. 

Figure 24 presents similar spectra and time history 
comparisons for the 65,000-ft.-altitude maneuver. In 
this case a 1-ft length of tube is used in conjunction 
with the 50 Hz Butte rworth filter. The output time de- 
lay, unlike the 8-ft tubing data (Fig. 20) is extremely 
small. Cross-correlation data presented in Fig. 25 indi- 
cate that for angle of attack the time delay is on the or- 
der of 4 msec as opposed to 1 00 msec for the 8-ft tubing 
sensor geometry. Thus, using a single pressure sensor 
geometry, coupled with proper low-pass filter signal 
conditioning, the HI-FADS system demonstrates qual- 
ity acoustical performance over a wide altitude range. 

Concluding Remarks 

The accuracy of a prototype nonintrusive airdata 
system designed for high-angle-of-attack measure- 
ments was demonstrated for quasi-steady maneuvers 
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up to 55° angle of attack during phase one of the F-18 
high alpha research vehicle (HARV) flight test pro- 
gram. The system evaluated in this paper consists of a 
matrix of nine pressure ports arranged in annular rings 
on the aircraft nose, and estimates the complete airdata 
set using flow modeling and nonlinear regression. The 
system performed well for steady or moderate maneu- 
vering flight conditions and computations were per- 
formed at a rate of 25 samples/sec. Satisfactory fre- 
quency response was verified up to the system Nyquist 
frequency of 12.5 Hz. For higher frequencies there is 
some question as to how well the high-angle-of-attack 
flush airdata sensing (HI-FADS) system would per- 
form. This paper addresses one of the primary reasons 
for this concern: the effects of acoustical distortions 
within the individual pressure sensors of the HI— FADS 
pressure matrix. 

To quantify these effects, a dynamic model which 
describes acoustical distortion was developed and 
solved in closed form for frequency response. The 
model was briefly verified by comparison to lab data. 
Simple design criteria which describe the dynamics of 
the primary harmonic were developed from the model. 
Use of these criteria allows for the geometry of the 
pressure measurement system to be tailored to a par- 
ticular frequency band or operating altitude. 

Using an input pressure data set, synthesized from 
HI-FADS and surface pressure measurements, the ef- 
fects of various internal sensor geometries and pre- 
sample filtering were analyzed and qualitatively eval- 
uated. Analyses and data presented indicate that the 
sensor frequency response characteristics vary greatly 
with altitude, thus it is difficult to select any single sen- 
sor geometry which gives acceptable acoustical perfor- 
mance at all altitudes. 

The altitude variability problem is overcome by us- 
ing presample signal conditioning to eliminate the res- 
onance effects. Caution must be used in scaling the 
transducer range to accommodate resonance. If the 
resonance causes the transducer to saturate, then when 
the transducer output is run through the low-pass fil- 
ter, the nonlinear distortion caused by the saturation 
process will result in a gain offset. Consequently the 
measured signal will be biased. Furthermore, cau- 
tion should be exercised if multiplexed or electroni- 
cally scanned pressure transducers are to be used in 
the nonintrusive airdata system. If multiplexed sen- 
sors such as the electronically scanned pressure (ESP) 


module are to be used, then the scan rate must be high 
enough to clearly capture any significant tubing reso- 
nances. Any harmonic which lies beyond the Nyquist 
frequency of the scan will be aliased and appear at a 
lower frequency within the band of interest 
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Fig. 2 Schematic of HI-FADS pressure sensing system 
internal geometry. 
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Fig. 3 Comparison of accoustical model frequency 
response to laboratory data for HI-FADS sensor ge- 
ometry, 2300-ft altitude. 
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Fig. 4 Schematic of spectral convolution scheme. 
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Fig. 5 Schematic of spectral deconvolution scheme. 
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Fig. 6 Calculated frequency response of typical Fig. 7 Calculated sensor frequency response; L = 1 ft, 

HI-FADS pressure sensor, L = 8 ft, D = 0.06 in., H p = D = 0.06 in., V = 0.01 in 3 , H p = 20,000 ft. 

20,000 ft. 
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Fig. 8 Calculated sensor frequency response; L = 8 ft, 
D = 0.06 in., V = 0.01 in 3 , H p = 65,000 ft 
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Fig. 9 Calculated sensor frequency response; L = 1 ft, 
D = 0.06 in., V = 0.01 in 3 , H p = 65,000 ft. 
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Fig. 10 Calculated sensor frequency response; L = 
1 ft, D = 0.06 in., V = 0.01 in 3 , H p = 20,000 ft; cou- 
pled with second-order, 50 Hz Butte rworth filter. 
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Fig. 11 Calculated sensor frequency response; L = 
1 ft, D = 0.06 in., V = 0.01 in 3 , H p = 65,000 ft; cou- 
pled with second-order, 50 Hz Butterworth filter. 
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Fig. 12 Complementary filtering scheme for merging Fig. 13 Time history of 20,000-ft altitude pushover- 
pressure data. pullup maneuver. 
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Fig. 14 Comparisons of input and output power spectra: ESP101 pressure and HI-FADS angle of attack. 
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Fig. 15 Time history comparisons of input and out- 
put HI-FADS parameter values, angle of attack, and 
Mach number. 
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Fig. 16 HI-FADS angle-of-attack time delay; corre- 
lation of input and output values. 
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(a) ESP101. (b) Angle of attack. 

Fig. 17 Comparisons of input and output power spectra: ESP101 pressure and HI— F/.DS angle of attack; L = 1 ft, 
H p = 20,000 ft. 
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Fig. 18 Time history comparisons on input and out- Fig. 19 Time history of 65,000-ft altitude simulated 

put HI-FADS parameter values: angle of attack, and pushover-pullup maneuver. 

Mach number, L = 1 ft; H p = 20,000 ft. 
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(a) Time history. (b) PSD. 

Fig. 20 Comparisons of input and output time history and power spectra: HI-FADS angle of attack; L = 8 ft, H v = 
65,000 ft 



Fig. 21 HI-FADS angle-of-attack time delay: correlation of input and output values; L = 8 ft, H p = 65,000 ft. 
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(a) ESP101. (b) Angle of attack. 

Fig. 22 Comparisons of input and output power spectra: ESP101 pressure and HI-FADS angle of attack; L = 1 ft, 
H p = 20,000 ft, with 50 Hz antialiasing filter. 
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Fig. 23 Time history comparisons of input and output (resonated) HI-FADS parameter values, angle of attack, 
and Mach number, L = 1 ft, H p = 20,000 ft, with 50 Hz antialiasing filter. 
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Fig. 24 Comparisons of input and output time history and power spectra: HI-FADS angle of attack; L = 1 ft, H p = 
65,000 ft, with 50 Hz antialiasing filter. 



Fig. 25 HI-FADS angle-of-attack sensor time delays: correlation of input and output values; L - 1 ft, H p 65,000 ft, 
with 50 Hz low-pass filter. 
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